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Abstract 


Dynamic  equations  are  developed  for  rubber  viscoelasticity  based 
upon  a  stick-slip  continuum  molecular  based  model.  The  model  con¬ 
sidered  is  a  continuum  simulation  of  a  tube  reptation  model  in  which 
a  chemically  cross-linked  (CC)  system  of  molecules  act  as  a  constraint 
box  per  unit  volume  for  a  physically  constrained  (PC)  system  of  mole¬ 
cules.  The  CC-system  carries  along  the  PC-system  during  instanta¬ 
neous  step  deformations.  The  subsequent  relaxation  of  the  PC-system 
is  determined  by  the  history  of  the  CC-system  and  thermodynamic 
considerations.  At  the  same  time  the  PC-system  deformation  acts  as 
an  internal  variable  affecting  the  relaxation  of  the  CC-system.  Mod¬ 
eling  of  this  relaxation  process  and  subsequent  employment  of  energy 
density  functions  constitutes  one  basis  to  model  viscoelastic  effects  in 
rubber  deformations. 


1  Introduction 


Various  molecular  and  phenomenological  models  have  been  developed  to 
model  both  small  and  large  deformations  in  rubber.  For  a  recent  review, 
see  the  survey  by  Johnson  [Jj.  To  provide  a  more  accurate  and  realistic 
model  of  the  deformation  and  relaxation  process  in  rubber  molecular  the¬ 
ory  [S],  [T]  and  continuum  mechanics  [0],  it  is  useful  to  integrate  various 
approaches. 

In  this  paper  we  use  the  molecular  models  of  Doi/Edwards  [DE]  and  John- 
son/Stacer  [.JS]  to  derive  a  class  of  nonlinear  distributed  parameter  systems 
(partial  differential  equations)  with  internal  strain  dynamics  that  include  the 
pseudo-phenomenological  models  of  [BPP],  [BPPGY],  These  latter  models 
have  provided  good  agreement  with  both  quasi-static  and  dynamic  data  for 
rubber  in  uniaxial  tension  and  in  shear. 

The  dynamic  model  for  rubber  viscoelasticity  presented  here  is  based  on 
the  continuum  simulation  [JS]  of  a  tube  reptation  model  considered  in  [DE], 
In  [DE]  step-strain  relaxation  of  polymers  is  modeled  with  constraint  (stick- 
slip)  theory  in  which  PC-molecules  deform  with  CC-molecules  during  a  large 
step-strain.  Then  the  PC-molecules  contract  and  creep  to  return  to  a  low 
energy  and  higher  entropy  state.  As  a  result,  the  total  energy  density  at  a 
constraint  strain  dissipates  in  time  and  a  viscoelastic  theory  results.  These 
models  are  based  on  cross-linking  rubber  network  theories  for  rubbers  and 
other  polymers  [S],  [T], 

A  continuum  molecular  model  of  rubber  viscoelasticity  is  proposed  in 
[JS]  and  entails  consideration  of  the  chemically  cross-linked  molecules  as 
providing  cells  or  boxes  with  entrapped  molecular  segments.  The  model 
involves  placing  a  unit  cell  or  box  at  each  point  of  the  rubber  continuum  and 
deriving  subsequent  equations  for  the  associated  principal  stretches. 


2  The  Continuum  Model  of  Johnson  and  Stacer 

In  this  section  and  the  next  we  summarize  the  arguments  of  Johnson  and 
Stacer  in  [.JS].  Let  the  CC  constraint  tube  have  length  L{t)  and  the  entrapped 
PC-molecule  have  length  £{t).  The  model  under  consideration  assumes  that 
the  length  of  the  CC  tube  is  approximated  as  a  step  function  of  time  and 
£(0)  =  I.  ( 0 ) .  Suppose  an  instantaneous  tensile  step-strain  deformation  at 
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time  of  the  CC-system  results  in  the  PC  molecule  having  length  £* ,  and  a 
model  of  the  PC-molecule  returning  to  its  original  contour  length  is  given  for 
t  >  0  by 

£{t)  =  1(0)  +  [t  -  £{0)]e-t/r. 

If  the  instantaneous  stretch  at  a  time  t0  for  the  CC-tube  is  A T0,  and  that  for 
the  PC  molecule  is  A 4 .  and  if  we  relate  the  stretches  in  the  form  A2-  =  I2-, 
we  hnd  the  relation 

Ai° =  (t^)  ALo ' 

Then,  in  the  time  interval  t0  <  t  <  ti  we  have 

£(t)=£0  +  ^AL0e~^^,  (1) 

To 

where  r  is  the  relaxation  time  for  the  slip  motion.  Next  we  suppose  that 
the  CC-box  is  subjected  to  further  instantaneous  tensile  step-deformations 
of  magnitudes  A  4  at  times  /.*.  /  —  1.2....  .  Setting  t  =  A,  we  hnd 

£(h)  =  4  +  ^A  T0e-(tl-to)/T, 

To 

Ti  =  To  +  ATo- 

For  t  in  the  time  interval  (4, 4),  we  have 

/(/)  =  4  •  [(/,  •  AC)  4K(! 

=  4+  -p-AT0e_(tl_t°l/T  -f  -^-ATa  e-h-T)A 
T0  Ti 

=  4  +  ^Af0e-(t-‘o)/T  +  -^Ah,e  ('  'l)/r.  (2) 

To  Ti 

Similarly  for  t2  <  t  <  t3  we  have 

1(4  =  4  +  ^AI0e_(1_to)/T  +  y-AT1e_(t_tl^T  +  ^Af2e_lt_i2)/T. 

Lo  L\  L 2 

Passing  to  the  limit  as  At*  =  ti  —  tt-_i  — »  0,  we  are  thus  led  to  the  formula 


<(<)  =  4  +  f  (W  dLjs>e-«-V’ds 

Jo  L  S  GCS 
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for  the  length  of  the  PC-molecule  at  any  time  t  in  terms  of  the  previous 
lengths  l(s),L(s)  and  the  rate  of  change  In  differential  form  this 

becomes 


d!  _  hi  / 1  1  dL\ 

dt~~~  \t~  L~dt  )  '  U 

3  Strain  Energy  Density 

To  use  the  stick-slip  model  in  continuum  simulation  of  the  reptation  model 
of  rubber  elasticity  one  considers  a  network  of  cells  or  boxes  in  the  rubber 
continuum  with  sides  Ai,  A2,  and  A3.  The  CC-box  will  have  positive  strain 
energy  density  Wcc(Aic,  A2c,  A3c)  for  all  stretches  except  when  Aic  =  A2c  = 
A3c  =  1.  Next,  a  box  for  the  PC-system  with  sides  parallel  to  those  of 
the  CC-box  is  defined  by  sides  Aip,  A2p,  A3p  along  with  an  energy  density 
fhpc^Aip,  A 2p,  A3p). 

Stresses  are  calculated  from  the  strain  energy  density  [VL]  by  determin¬ 
ing  how  the  energy  density  function  changes  with  respect  to  changes  in  the 
stretches  or  displacements.  In  the  model  considered,  the  strain  energy  density 
of  the  rubber  continuum  is  assumed  to  have  the  form 


W  =  hPcc  +  ITpc 

=  W4  A  I.  :  A2c,  A;v  I  +  IPpc(  Alp,  A 2p,  ^3  p)-  (4) 

To  hnd  the  stresses  at  a  generic  point  of  the  rubber  continuum  one  must 
determine  how  W  changes  with  respect  to  stretches  of  the  CC-system.  The 
stretches  of  the  PC-boxes  are  then  treated  as  internal  variables  depending 
on  the  stretches  of  the  CC-system. 

The  Cauchy  stress  in  the  principal  direction  e?,  where  ej  is  a  unit  vector 
in  the  Xj  direction,  is  given  by 

dW 

T ’  =  ^ax~  -  p  (5) 

where  P  is  the  hydrostatic  stress.  If  we  consider  (4)  with  the  Ajp  as  internal 
variables  depending  on  the  A 4C,  (5)  becomes 
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_  .  d\Vcc  *dWpcd\lp 

Tj  —  -W  si\  T  A'c  *• 


dX 


JC 


i— i  $ X{p  OXjc 


(6) 


In  obtaining  (3)  we  have  related  the  instantaneous  tensile  step  deforma¬ 
tions  Aik,  A Lj.  at  t  —  4  according  to  the  formula 

A4  _  4 

AL  k  L  k 

Using  this  observation  in  the  analogy  of  A4,A4*,  with  principal  stretches 
XiPi  A ic,  respectively,  we  may  write 


^ jp  c 

dXir  A  ir  jt' 


(7) 


In  the  case  of  tensile  deformations,  choosing  j  =  1  for  the  direction  of 
loading  and  using  (7)  and  (4)  in  (6),  we  have 


n  ,  dW*  ,  ,  dWm  n 

t~2  —  t3  —  0  —  A2c  — - h  X2p— - r. 


d\ 


2c 


A 


2  v 


(8) 


With  the  hydrostatic  stress  determined  from  (8),  the  tensile  Cauchy  stress 
is  thus  given  by 


Tl 


(.  dwcc  dwcc\  (  dwpc  dwpc\ 

\  lc  dXlc  2c  dX2c  )  V  1P  d\i P  2P  dX2p  )  ' 


(9) 


4  A  Dynamic  Continuum  Model 

To  demonstrate  how  the  box  model  detailed  above  can  be  applied  to  the 
analysis  of  rubber  undergoing  large  dynamic  tensile  strains,  we  consider  the 
strain  energy  function  suggested  by  Johnson  and  Stacer  (based  on  experi¬ 
mental  data  of  Young  and  Danik  -  see  [JQYD])  given  by 


Wc4  =  105.0(/i  -  3)  +  103. 0(4  -  3), 
Wpc  =  169.0(/i  -  3)  +  0.0138(4  -  3)2  +  7.89(4  -  3)3, 

where  4  =  Xjc  +  A2C  +  A2C,  and  4  =  A2CA2C  +  XPcX^c  +  A2CA2C. 


(10) 
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We  also  impose  the  incompressibility  condition  AiA2A3  =  1  on  the  prin¬ 
cipal  stretches  A,.,  =  (A A2c,  A3c)  and  \pc  =  (A|;t..  A2p,  A3p). 

From  (9)  the  engineering  stress  ny  =  ti/Aic  is  given  by 


_  dWcc  A2c  dwcc  |  Alp  dwpc  X2p  dWpc 

d\\c  Aic  d\ 2c  Aic  dX\p  Aic  dX2P 

where,  from  (10),  we  have 


(11) 


dwcc  a2c  dwcc 

dXj c  Aic  dX2c 

=  210Alc  +  206(AlcA2c  +  AicA2c) 

—  - — (210A2c  +  206(A2cA2c  +  A2cA2c)) 

Me 

=  210Aie  +  206Aic A^c  +  20()A),  A3,  - 
210AL  206ALAL  206A|oA2c 


For  the  tension  in  the  incompressible  CC-system  we  see  that 


A 


2c 


—  A3c  — 


1 

y/Xic 


Then,  from  (12),  we  can  compute 


dwcc  X2c  dwcc 

dXi c  Aic  <9A2c 

210Alc  +  206 


210 

AL 


206 

AfT' 


Next  we  deal  with  the  PC-molecules  and  find 


(12) 


(13) 


\Vpc  -  169.0  (A,p  +  \\p  +  \lp  3]  0.0318  (A^p  +  \\p  +  \\p  3]  I 


7.89  {\\pX\p  +  AjpAgp  +  A^A 

Using  the  relationship 

X2p  —  A3p  = 


2  -  2 

3  p  ° 


(14) 
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for  tension,  we  obtain  from  (14) 


Aip  d\ bpC  A2p  d\ bpC 

Ale  ^Aip  Ale  dX2p 

338.0  +  0.0552  (  2Xlp  +  -L 

\  Aip 


(15) 


Letting  uc  represent  the  deformation  of  the  CC-box  we  have 


Ale  1  ~\~ 

A“2  Ri  1  -  2dxiic 

A^c3  ~  1  -  3 dxuc 


while 

Then 


Aip  1  +  dxiip 


dWcc 

dXic 


A 2c  dWcc 

Ale  dX2c 


fa  210(1  +  dxuc)  +  206  -  210(1 

-206(1  -  3dxuc) 

=  1248 dxuc. 


2dxuc) 


(16) 

(17) 


Let  up  represent  the  deformation  of  the  PC-box.  Then,  proceeding  as 
above 


Alp  C^Lfpe 
Ale  dX\p 


X2p  dWm 
Aip  dX2p 


1014(1  -  dxuc)dxup. 


Thus,  using  (11)  we  obtain  the  engineering  stress  approximation 


(18) 


u  fa  1248 dxuc  +  1014(1  —  dxuc)dxup 
=  (1248  —  1014 dxUp)dxuc  +  10143^^. 

In  (3)  taking  Ai p(t)  for  £(t),  Alf:(/,)  for  L(t)  and  using  l0  =  1(0)  =  Aip(0) 
we  find 


(19) 
=  1, 
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(20) 


d\\p  1 

dt  t 

Using  the  definitions  Alp  =  1  +  dxup ,  A1(:  =  1  +  dxuc  in  (20)  we  have 


1  1  d\lc 

t  Aic  dt 


3  f  32v  1 

—  (1  +  dxup)  =  1/t  -  <  1/t  -  (1  +  (Irti.j  '  — j  (1  +  (21) 

In  terms  of  the  infintesimal  strains  cj  =  c^itp  for  the  PC  system  and  £  =  dxuc 
for  the  CC  system,  we  have 

T£  i  +  ci  =  Tt(l  +  £i)/(l  +  £)•  (22) 

If  we  use  the  approximation  (1+t)-1  ~  1  — £  +  .  .  (22)  provides  a  relationship 

for  “internal  dynamics” 


T£ l  +  Cl  =  T£(l  +  £i)(l  -  c). 


(23) 


Ignoring  the  terms  ££!£  in  (23),  we  obtain  a  further  approximation  of  the 
internal  dynamics  given  by 


£i  +  —£\  =  £  +  £  [ci  -  e] . 


(24) 


We  remark  that  in  integrated  form  the  equation  (24)  can  be  written  (using 
an  integration  by  parts  after  integration) 


£i(t)  =  s(t) — £(0)e  4 - f  e(s)e  d  '  ds-\-  [  [%(a)  — £(a)]  e(a)e  ^  S^’ds. 

T  JO  JO 

(25) 

In  terms  of  the  deformations  3xup,  3xuc  this  is  the  same  as 


3xv,p 


3xuc - f  3xuc[s,x)e  d  ’  ds 

T  Jo 

,  [\-A  <  5  a  /  ^d2uc(s,x)  ,,  s]/T  , 

+  J^[OxUp{s,x)  -  axUc(s,X) \ — — e  ds 

-3xuc(0,  x)e~tlT.  (26) 
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We  next  consider  the  longitudinal  vibration  of  a  rubber  rod  with  (unde¬ 
formed)  cross-sectional  area  Ac  and  mass  density  p.  Using  (19)  we  write  the 
equation  of  motion  (recall  that  S  =  Acu  is  the  engineering  stress  resultant  - 
see  (2.1)  of  [BPPGY]) 


pAJ)fv,  ~  dx[Ac(  1248  -  lQUdxup)dxuc\  -  dx{10UAcdxu^  =  q,  (27) 

where  q  is  the  applied  load  (in  force/unit  length)  and  this  equation  must  be 
solved  with  (26)  (or  equivalently  (22)),  or  one  of  its  approximations  (23)  or 
(24)  to  provide  a  defining  coupling  for  dxup  and  dxuc. 

In  summary,  the  full  nonlinear  model  we  have  derived  for  the  CC-PC 
molecular  system  in  tensile  deformations  for  a  rod  have  the  form 

pdfuc  -  dxa(6,6\)  =  q/Ac  (28) 

<t(c,  £,)  =  (1248  -  lOPlcOc  +  1014c,  (29) 

T 2  |  '  £j  —  T£  (  1  ’  c  |  ) /  ( 1  •  c )  (30) 

with  £  =  dxuaiG\  =  dxup  and  q  the  applied  load  (in  linear  force  density). 

We  remark  that  the  system  (27),  (23)  can  be  viewed  as  a  member  of  a 
class  of  systems  similar  to  those  treated  in  [ABP]  and  the  methods  presented 
there  can  be  used  to  guarantee  well-posedness  (existence  and  uniqueness)  of 
the  corresponding  variational  forms  of  this  system.  The  models  treated  in 
[ABP]  can  be  properly  viewed  as  generalizations  of  (27),  (23)  that  permit 
nonlinear  instantaneous  elastic  responses  as  well  as  other  refinements. 


5  Qualitative  Analysis 

In  this  section  we  turn  to  a  brief  discussion  of  some  of  the  qualitative  prop¬ 
erties  of  the  model  derived  in  the  previous  section.  We  first  linearize  (29)  to 
obtain 


<rL(e,e)  =  1248s  +  1014£l.  (31) 

We  use  the  linearized  version  of  (30)  -  essentially  (24)  in  which  we  ignore  the 
term  £(£i  —  e)  -  given  by 
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(32) 


1 

£l  ~\ - £  |  —  £. 

T 

Integrating  and  using  e(0)  =  5a.wc(0,  x )  =  0,  we  obtain 


£i(t)  =  £{t)  —  —  f  £(s)e  ^  s^Tds.  (33) 

r  Jo 

This  results  in  the  linearized  system 

pd\uc  —  2262 dxuc  -\ — - — -  [  dluc(s ,  -)e~^~s^Tds  =  q/Ac.  (34) 

T  Jo 

Making  the  change  of  dependent  variable  given  by 


we  obtain 


2262 

v  =  - uc 

P 


1014 

pr 


1014  1248  2262 

Vff  H - Vf - -v - vxx 

2262 r  2262t2  p 

.  1248  ,n  ft  1248  '  1  . 

+W  l 


(35) 


(36) 


where 


2262 

S=JA?- 


(37) 


Next  consider  the  standard  Sturm-Liouville  problem  resulting  from  sep¬ 
aration  of  variables  in  systems  such  as  equation  (36).  The  generic  problem 
is  given  by 


-y"  ~  Ay  =  0 

y(a)-h0y'(a )  =  0,  h0  >  0 

y(b)  +  hiy'(b)  =  0,  hi  >  0. 


(38) 

(39) 


For  this  problem  it  is  known  [CL]  that  there  is  a  sequence  of  eigenvalues 
0  <  Ai  <  A2  .  .  .  <  X-n  <  ■  ■  ■  /*  00  and  corresponding  complete  family  of 
orthonormal  eigenfunctions  (pi,  (^2,  •  •  ■ ,  ‘pn,  ■  ■  ■■ 
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(40) 


From  this  one  is  led  to  consider  a  solution  of  (36)  in  the  form 


OO 

v{t,x)  =  Vn(t)ipn(x). 

n= 1 


Also  expanding  the  forcing  function 


/(M)  =  fn{t)<Pn{x), 

n= 1 


we  formally  obtain  from  (36)  the  equations 


1014 


W  +  + 


+ 


2262 r 
/  1248  \2  A 


/2262An  1248  \ 
V  -p  2262r2  J 

1 


■(f) 


j  jo  e  s)(»„(s)  +  -un(s))ds  =  /„(*). 


Letting 


we  obtain, 


\2262r 

(  1248  \2  A 


=  (22627)  /o  + 


/  1248  \ 


1248 


=  I22627J  +  “  2262t W"(^' 


(41) 

(42) 


(43) 


We  write  (41)  as  first  order  vector  systems  by  defining  new  variables 


«t) 

«;(() 

«3<0 

4 


A 


n 


®n(^) 5 
(^)  5 

(^) : 

1248 

2262r 

2262  x  1248 
«  ”  ~~  2262r2 ' 


Then,  from  (41)  we  have  the  equivalent  systems 


d 

dt 


/ «?  ^ 

/ 

= 

V  ; 

V 

0 
■  A 

d2 


1 

-1014 

2262t 

d2 


0 


-1 

-d 


/  <  \ 

( 

0  ^ 

«2  + 

fn 

V  <4  I 

l 

0  ) 

(44) 


(45) 
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Letting 


An  =  -A 


1  0 

-1014 


-1 


n  2262 r 

&  d2  -d 


(46) 


we 


find 


\ci-An\  =  e+(—+d)e  +  (—d+<p+A1)ic 

ls  n|  s  V2262 r  7  S  \2262r  7  S 


(  1014 


A  4  ^ 

+A nd  -\ - . 

T 


(47) 


The  Routh-Hurwitz  Theorem  [G]  tells  us  that  the  roots  of  the  polynomial 
equation 


23  +  a\z 2  +  a^z  +  «3  =  0 


have  negative  real  parts  if  the  following  conditions  are  met 

(i)  «i  >  0 


(ii) 


a  |  a3 
1  «2 


>  0 


(iii)  a3  >  0. 

Thus,  the  roots  of  |(  /  —  /I n  |  =  0  have  negative  real  parts  if 


1014 

2262t 


d 


1 


1014 

2262t 


A  nd  +  4" 
d  +  d2+Ar 


>  0. 


Thus,  the  eigenvalues  of  A „  have  negative  real  parts  if 


d 


Now,  recalling  that 


(  1014  A2  1014  d2  n 

+  d  +^^An - >0. 


,2262r 


2262r 


d  = 


A  n  — 


1248 

2262r 

2262A, 

P 


1248 
2262r2 ’ 


(48) 


(49) 
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we  find  that  simple  calculations  reveal  that  inequality  (49)  is  equivalent  to 


\n  >  0. 


(60) 


Since  condition  (50)  holds  for  the  eigenvalues  of  (38), (39),  we  see  that  the 
eigenvalues  of  An  have  negative  real  part  and  hence,  the  zero  solution  of 
(36)  with  /  =  0  is  asymptotically  stable.  That  is,  all  modes  of  (36)  are 
exponentially  damped. 

From  (42),  (44),  (45),  and  (46)  we  have  that  the  temporal  coefficients  of 
(40)  satisfy 


/ 

Vn{t)  \ 

( 

MO)  ^ 

1  /•  t  , , 

( 

°  ^ 

Vn{t) 

--  etAn 

l’n(O) 

I  T  / 

fn{s) 

\ 

wn(t )  ) 

0  j 

J  Jo 

0  J 

where  for  each  n, 


1  M°) 
M°) 

\  0 


oo. 


(51) 


6  Vibration  of  a  Rubber  Rod  with  a  Tip  Mass 


In  this  section  we  consider  the  Lagrange  formulation  for  longitudinal  vibra¬ 
tions  of  a  rubber  rod  fixed  at  x  =  0  with  a  tip  mass  at  x  =  b.  For  this  we 
return  to  the  linearized  equation  (34),  set  q  =  0  and  adjoin  the  boundary 
condition 

d2u 

M^(t,  b )  =  FU  +  F(t)  +  Mg,  (52) 

where  M  is  the  tip  mass,  g  the  gravitational  constant,  F  is  the  force  applied 
to  the  tip  mass,  and  S  —  Ac(Tl  is  the  internal  (engineering)  linearized  stress 
resultant  given  by,  in  this  case, 

S  -  2262 Acdxuc  -  10UAc  f  dxuc{Sl  x)e~^/Tds.  (53) 

r  Jo 


Using  (35)  we  obtain 


S  -  pAcdxv , 


(64) 
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and  from  (34) 

d^Uc  =  vxx  -\ - —  q.  (55) 

pAc 

For  the  problem  of  the  vibration  of  a  rod  with  a  tip  mass,  as  indicated  in  the 
first  paragraph  above,  we  have  q  =  0.  Thus  from  (55)  we  have 

d]uc  =  vxx.  (56) 

Moreover  from  (52),  (54)  and  (56)  we  have 

Mvxx(t,  b)  =  - pAcvx(t ,  b )  +  F(t)  +  Mg.  (57) 


Formally  employing  the  Sturm- Liouville  theory  and  series  solution  from  the 
previous  section  in  (57),  we  have 


E 


M\t,if„(b) 


-  pAc 


vn(t)  =  F(t.)  +  Mg. 


(68) 


We  may  then  formally  write 


OO 


F(t)  +  Mg  =  +  Mg,<pn)ipn, 

n= 1 

(69) 

where  (,)  is  the  L2— inner- product.  Then,  using  (58)  and  (59),  and  the  fact 
that,  under  appropriate  regularity  assumptions,  the  series  in  (59)  converges 
pointwise  for  every  x  in  (0,  6),  we  have 

(~M\n  -  YpAe)  Vn (f)  =  +  Vfl-.r") 

(60) 

or 

(F{t)  +  Mg^n) 

”"(i)  =  -MK  -  &A.  ■ 

(61) 

Thus, 

u  y  {F{t)  +  Mg,ipn) 

(  -M\n-±pAcVn(X>' 

(62) 

We  have  used  formal  arguments  with  (57)  and  (58),  enabling  us  to  write 
(61)  for  vn(t),  and  thus  (62)  for  v(t ,  x ).  However,  having  this  form  for  v(t,  x ), 
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we  can  invert  our  arguments  to  give  a  rigorous  development  for  this  rep¬ 
resentation  (62).  Specifically,  using  known  facts  about  completeness  of  the 
eigenfunctions  of  Sturm- Liouville  problems  and  convergence  properties  of  ex¬ 
pansions  in  terms  of  the®  eigenfunctions  we  can  establish  (62)  with  complete 
rigor.  Indeed,  taking  a  =  0,  h0  =  0  in  (38),  we  can  readily  compute  the 
eigenvalues  and  eigenfunctions  and  validate  our  formal  steps. 

Returning  to  (35)  we  can  solve  for  u,  giving  the  displacement  of  the 
rubber  rod.  Thus, 


u, 


.(t, 


X  = 


2262 


[f  1248  / j  )  (  1  \  1248  j  . 

J  e  2262  V  vl  Vg  -| — vj  as  +  e  lUc(U,  x) 


P  MS-t  ,/n  ^  ,  1014  p 


2262 

-1248 


e  2262x^(0,  X)  + 


2262  2262t 


f  e  2262r  0  S)V(s:x)ds 

Jo 


+e  2262r  tuc(0,  x). 


Thus,  using  (62)  we  have 


u 


.{t, 


x 


P 


+ 


2262 

1014 


-  1248  t  <^(0)  +  Mg ,  <pn)  . 
e  2262T  - rr1^^) 


i=i  M\n  pAc 


P 


2262  2262r 


['  Y  (r (.«)  +  Mg,  y.)  .  w 

Jo  h  -MK  - 


_  1248  ,  , 

+  e  22621  Wc(0,  X). 


(63) 


We  see  that  (63)  is  an  explicit  formula  for  the  displacement  of  the  rubber 
rod.  If 

\F(s)\  eA^Fs  <  Ce~ts ,  e  >  0, 

we  see  from  (63)  that,  asymptotically,  only  the  weight  of  the  tip  mass  deter¬ 
mines  the  displacement  of  the  rod. 


7  Conclusion 

We  have  presented  a  dynamic  simulation  of  rubber  viscoelasticity  based  on 
reptation  models.  Deformed  molecular  chains  or  segments  entrapped  between 
cross-linked  molecules  or  molecular  chains  tend  to  return  to  their  positions 
prior  to  deformation.  This  is  due  to  the  fact  that  their  original  positions 
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are  positions  of  lower  energy  and  higher  entropy.  However,  due  to  the  phys¬ 
ical  and  electrostatic  barriers  created  as  a  result  of  new  configurations  and 
conformations  the  entrapped  and  strained  molecules  only  creep  to  their  orig¬ 
inal  confirmation.  The  model  presented  adheres  to  these  observations  and 
also  points  out  the  interaction  of  density  and  relaxation  time  on  the  creep 
behaviour.  These  models  can  be  related  in  a  direct  manner  to  previously 
derived  (pseudo-phenomenological)  models  based  on  data,  from  quasi-static 
and  dynamic  experiments  with  rubber  rods  in  uniaxial  tension  and  in  shear. 
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